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2.1 LiDAR RANGING PRINCIPLE 


Ranging is a very important application of LiDAR systems. There are two main 
LiDAR ranging modes: pulse-based ranging and phase-based ranging. The ranging 
principles of the two modes are described in the following sections. 


2.1.1. RANGING PRINCIPLE OF PULSED LIDAR 


At present, most LiDAR systems work in the pulsed-based ranging mode. That is, 
the laser emits a very narrow pulse (the pulse width is usually less than 50 ns) toward 
the target, and then part of the pulse is reflected by the target and returns to the 
receiver. The LiDAR system records the time of flight (TOF) between the emission 
and reception of a laser pulse and then calculates the distance between the senor and 
the target according to the speed of light, as in Equation (2.1). Figure 2.1 shows the 
ranging principle of pulsed LiDAR. 


1 
R=—-c't 2.1 
5 (2.1) 
where R is the distance between the sensor and target, c is the speed of light in air, 
and ¢ is the TOF of the laser between emission and reception. 
Differentiating Equation (2.1), we have Equation (2.2): 


AR= 0° AN (2.2) 


where AR is the ranging resolution, representing the minimum distance at which two 
objects can be distinguished, depending on the time sampling interval. 

Additionally, in order to ensure that the returns of multiple laser beams can be 
distinguished, the pulsed LiDAR system must transmit the next laser pulse after 
receiving the return of the previous laser pulse. In this case, the maximum mea- 
surement distance (Rmax) needs to be considered (Lai, 2010). It is calculated by the 
longest time measurement (tnax), aS in Equation (2.3): 


Rmax = — C+ bax (2.3) 


In practical applications, the maximum measurement distance is affected by 
many factors, such as laser power, beam divergence, atmospheric transmission, 
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FIGURE 2.1 Ranging principle of pulsed LIDAR. 


target reflection characteristics, detector sensitivity, flight altitude, and deviation of 
flight attitude record. 


2.1.2 RANGING PRINCIPLE OF PHASED LIDAR 


The phase-based LiDAR is also called continuous wave LiDAR. It transmits an 
intensity-modulated continuous laser which is reflected by the target and then 
returns to the receiver. The distance between the target and sensor is calculated by 
the phase variation of the laser during the round trip, as in Figure 2.2. Generally, 
the ranging precision of a phased LiDAR system is higher than that of a pulsed 
LiDAR. 

Assuming that the modulated continuous laser wave is a sine waveform with a 
period of T, the ratio of the time interval of the measured phase difference to the 
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FIGURE 2.2 Ranging principle of phased LIDAR. 


LiDAR Remote Sensing Principles 5 


period T is exactly the ratio of the phase difference of the transmitted wave and the 
received wave to the phase 27, as in Equation (2.4). 


t_ 9 

ee, 24 
T 20 Pe 
where @ is the phase difference of the transmitted wave and the received wave. 


According to Equation (2.4), the distance between the sensor and target (R) is 
expressed as in Equation (2.5). 


c T 
acancery 2.5 
2 2n ọ (2.5) 
where the period T is the reciprocal of the modulated frequency f of the sine wave. 
In this case, R can be expressed as in Equation (2.6). 


R=—-—.6 (2.6) 


where the ratio of the speed of light (c) to the modulated frequency (f) is exactly the 
wavelength (A). Hence, the distance R is as in Equation (2.7). 


nr 
R=7 rt (27) 


Differentiating Equation (2.7), we get the ranging resolution (AR) of phased 
LiDAR, as in Equation (2.8). 


AR= 7-9 (2.8) 


Equation (2.8) indicates that the ranging resolution of the phased LiDAR is 
related to not only the phase resolution but also the laser wavelength. In contrast, 
the ranging resolution of the pulsed LiDAR is only related to the time resolution 
[Equation (2.2)]. With the same phase (time) resolution, the ranging resolution 
of the phased LiDAR system could be improved by using a shorter wavelength. 
Hence, the phased LiDAR is usually used in applications with high ranging preci- 
sion. For example, for a continuous laser wave with the period of 1 second, assum- 
ing that the phase resolution is 0.4°, its ranging resolution is 0.1 mm. However, to 
achieve the same ranging resolution, the pulsed LiDAR system must have a time 
resolution of 1 picosecond, which requires higher standards for the time interval 
counting hardware system. 

Although the phased LiDAR system is more suitable for high-precision ranging 
applications, it is usually used for short-range detection. This is because it requires 
high power to emit a continuous laser. Also, it may cause ranging ambiguity. The 
maximum measured distance (Rmax) Of the phase-based LIDAR system which can 


max. 
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uniquely determine the target is the distance measured when the continuous wave 
travels exactly one period, as in Equation (2.9). 


Rmax = ` nax = N= (2.9) 
T 


Equation (2.8) indicates that the ranging resolution (AR) of the phased LIDAR 
system is related to the phase resolution (A@) and the laser wavelength (A). The 
shorter the laser wavelength used, the higher the ranging resolution. Hence, phased 
LiDAR makes it relatively easy to achieve high ranging resolution. However, 
Equation (2.9) shows that the maximum measured distance (R nax) of phased LIDAR 
is determined by the wavelength (À). The longer the wavelength, the longer the maxi- 
mum measured distance. Long-distance measurement and high ranging resolution 
both cannot be achieved in a continuous wave LiDAR system. Therefore, in order to 
meet the requirements of practical applications, the phased LiDAR system is usually 
equipped with a signal frequency modulation device that can adjust multiple fre- 
quencies. Additionally, high-power phased lasers are difficult to produce. Currently, 
most LiDAR systems use a pulsed ranging system. The LiDAR mentioned in the rest 
of this book refers to the pulsed LiDAR. 


2.1.3 RANGING PRECISION OF LIDAR 


Ranging precision is a pivotal parameter of a LIDAR system and is different from 
ranging resolution. Specifically, the ranging resolution refers to the limit of the 
ability to measure the distance between two objects along the same line of sight, 
while the ranging accuracy is the standard deviation of the target distance estima- 
tion in the presence of noise. The ranging precision is related to the signal-to-noise 
ratio (SNR) of the returned laser signal and the emitted laser pulse (Wehr & Lohr, 
1999). 

Usually, the ranging precisions of pulsed LiDAR and phased LiDAR are expressed 
as in Equations (2.10) and (2.11), respectively. 


o ho i (2.10) 
Rpulse 2 rise Phe . 
Boy 
ye (2.11) 
4a Pr, 


where O,,,,. IS the ranging precision of pulsed LiDAR, œ indicates a proportional 
relationship, fse is the pulse rise time, Bpuise is the noise width of pulsed LIDAR, Pr sai 
is the peak power of the received signal of pulsed LIDAR, Og.,, is the ranging preci- 
sion of phased LiDAR, / is the continuous laser wavelength, B,,, is the noise width 
of phased LIDAR, and Pp, is the average power of the received signal of phased 
LiDAR. 

Here we compare the performance of the two ranging modes with the assumption 
of measuring the same target. In this case, the received power can be replaced by 


LiDAR Remote Sensing Principles 7 


the transmitted power. That is, the peak power Pr of the received laser of pulsed 
LiDAR is replaced by the peak power of the transmitted pulse Prea» and the average 
power of the received laser of phased LIDAR Pr, is replaced by that of the transmit- 
ted laser P,- The ratio of ranging precisions of pulsed LIDAR to phased LiDAR is 
expressed as in Equation (2.12). 


O Ryuise Pry . Bpuise 


c 
cc 21> —- hise” 2.12 
(oy Rew A P, Tpeak Bew i i 
Assuming that: 
trise ° ; (2.13) 
Bogs 
Equation (2.12) can be rearranged as: 
o i 
Rpulse ee on P fe frise N Prw 214 
O Rw Boy Prex 


With the assumption of the pulse rise time of pulsed LIDAR of 1 ns, the peak power 
of transmitted pulse of 2000 W, the frequency of phased LIDAR of 10 MHz, the noise 
width of phased LIDAR of 7 kHz, and the average transmit power of 1 W, the ratio of 
ranging precisions of pulsed LiDAR to phased LiDAR can be calculated by Equation 
(2.14), that is about 0.0012. This indicates that although the transmit power of the 
pulsed LIDAR is 2000 times that of the phased LIDAR, the ranging precision of the 
pulsed LiDAR can theoretically only reach 85 times that of the phased LiDAR. 


2.2 LiDAR RADIATION PRINCIPLE 


The LiDAR system performs non-contact detection of an object’s surface by transmit- 
ting laser pulses and receiving the returned energy. Its remote sensing (RS) process 
includes four steps (Figure 2.3): (1) transmission of the emitted laser toward the target, 
that is, the pulse downward; (2) interaction of pulse and target surface, including inter- 
ception, scattering, absorption, etc.; (3) transmission of the scattered laser toward the 
sensor, that is, pulse upward; and (4) reception of the scattered laser by the receiver. 
Through this process, the received LiDAR signal is a comprehensive expression of 
the sensor configuration, the target surface characteristics, and the surrounding envi- 
ronment. In order to interpret the LIDAR signal, revert the LIDAR RS process, and 
retrieve the target information, it is necessary to understand the radiation principle of 
LiDAR RS. This section presents the radiation principle of the LiDAR system, which 
is also the core of LiDAR’s three-dimensional (3D) imaging mechanism. 


2.2.1 LIDAR Equation 


The physical process of the laser beam penetrating the atmosphere, interacting with 
scatterers, and returning to the receiver is very complex. Usually, beginners would 
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(4) Laser power received by the sensor: 


(1) Laser power emitted by the sensor: 


(3) Laser power scattered by the target: 


| ee 


FIGURE 2.3 LiDAR RS process. (P; Laser power emitted by the sensor; A,: effective area 
of target [scatterer]; A,,,.,: laser footprint area; p: reflectance of target [scatterer]; Q: scattered 
solid angle of target toward the receiver; R: distance between sensor and target; D,: diameter 
of LiDAR receiver; Nam: effect of atmospheric transfer; ys: efficiency of sensor transmitting, 


receiving, and processing signals.) 


simplify this complex process of LiDAR RS. The LiDAR equation is a physically 
based formula to compute the power returned to a receiver for given transmitted laser 
power, optical properties of the medium through which the LiDAR beam passes, 
and target properties. It is deducted based on two strong assumptions: (1) the power 
emitted by the laser is uniformly distributed within the footprint and (2) the scatterer 
uniformly disperses the incident power into a conical solid angle of Q (Wagner et al., 
2006). Based on the assumptions, the derivation process of the LiDAR equation is 
as follows. 
The footprint area of the laser beam at the distance of the scatterer is: 


TR’ B; 


Arse = 
4 


(2.15) 


where R is the distance between the sensor and scatterer and p, is the tan- 
gent of the laser beam divergence angle. Based on assumption (1) earlier, the 
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incident laser power density of the scatterer surface Ss is expressed as in 
Equation (2.16): 


4P, 
rR’ B? 


Ss = (2.16) 


where P, is the total emitted laser power. Equation (2.16) describes the property by 
which the laser power density gradually decays with the transmitted distance R. 
Considering the reflectance p and the effective area A, of the scatterer, the laser 
power dispersed by the scatterer P, is expressed as in Equation (2.17). 


4P, 
R? B? 


Based on assumption (2) earlier, if the scattered conical solid angle of the scat- 
terer overlaps with the field of view of the LiDAR receiver, the returned laser power 
density is expressed as in Equation (2.18). 


P= pA, (2.17) 


4P. 1 
mR? B? PA OR? 


S, = 


(2.18) 


Finally, considering the receiver aperture and the influence of the atmosphere 
and hardware systems, the LIDAR equation of a single scatterer can be expressed by 
Equation (2.19) as: 


a nD? 
7 PA oR $ 4 Natm "sys (2.19) 


Prit) = 

R, mae 
where Pp ;(t) represents the returned laser power from the i-th scatterer at the time t, 
D, is the diameter of the LiDAR receiver, Q is the scattered solid angle of the scat- 
terer, Nam is the effect of atmosphere transfer on laser power, and Nys is the effect of 
the sensor transmitting, receiving, and processing signals. 


All the parameters related to the scatterer can be combined into one, which is 
called the laser radar cross section o (LRCS) (Dai, 2002): 


ed g PA (2.20) 


Equation (2.20) shows that LRCS is determined by the effective area A,, reflec- 
tance p, and scattered solid angle Q of the scatterer. 

In the case where multiple scatterers are distributed in the propagated path of 
the laser, the received signal can be expressed as a general expression of the LIDAR 
equation, as in Equation (2.21). 


N 2 
D; 
Pum (t) = S tate 4nR# B? P, z Oi (R;) (2.21) 
i=l ree 
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where Pum (t) is the accumulation of the return signals of N scatterers at time t and 
the relationship of the time f and the distance R; is t= *~“, * represents the convolu- 
tion operation. 


2.2.2 LIDAR Waverorm MobEL 


Once there are multiple scatterers at different distances from the sensor in the laser 
propagation direction, the emitted laser might be reflected multiple times, and the 
reflected laser intensity and the traveled distance are both different for each time. 
The LiDAR waveform records the laser intensity according to the returned time 
delay, that is, laser intensity as a function of time. As in Figure 2.4, the LIDAR 
waveform can be quantified as a one-dimensional waveform with time as the hori- 
zontal axis and intensity as the vertical axis. It is theoretically feasible to inverse 
the characteristics of scatterers by accurately interpreting the LiDAR waveform. 
However, many parameters in the LIDAR equation are difficult to obtain accurately, 
such as Nam and MNsys- Also, the parameters such as Q and p are different for dif- 
ferent objects and under different observation conditions. This makes it very dif- 
ficult to directly inverse the surface characteristics from the LiDAR waveform. 
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FIGURE 2.4 Schematic diagram of LIDAR waveform. 
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Therefore, in practical applications, the LiDAR equation need to be further simpli- 
fied based on some specific assumptions. 

The LiDAR waveform is a product of the interaction between the laser pulse and 
Earth’s surface. The laser-surface interaction process can be expressed as the trans- 
fer function. Since the emitted laser pulse is close to a one-dimensional Gaussian 
distribution in the time domain, the LiDAR waveform is just the convolution of the 
emitted pulse and the transfer function in the time domain, which is usually approxi- 
mated by a mixed Gaussian model. That is, the LIDAR waveform is regarded as a 
superposition of multiple Gaussian functions, as in Equations (2.22) and (2.23). 


N 
Pan = > fi) (2.22) 
i=1 


(t-u)? 


fiQ=Ae 2° (2.23) 


where f;(t) is the i-th Gaussian component and A;, U;, and O; are the amplitude, 
mean, and standard deviation of the i-th Gaussian component, respectively. Based 
on this principle, the mixed Gaussian function can be used to decompose the wave- 
form, and the characteristic parameters of the ground object (such as reflectivity 
and geometric shape) can be obtained by inversion. For more details on waveform 
decomposition, see Section 4.2 of this book. 


2.2.3 LiDAR RapDIATIVE TRANSFER MODEL 


In addition to the LiDAR equation and the LiDAR waveform model based on very 
strong assumptions, some researchers have proposed relatively complex LiDAR 
radiative transfer models (RTMs). The RS radiative transfer models are considered 
celebrated tools for modeling, validating, and interpreting the RS signals, as they can 
simulate the interaction mechanisms of electromagnetic waves and Earth’s surface 
that give rise to the signals. 

To date, many researchers have developed a wide range of LIDAR RTMs with 
diverse complexity and efficiency. According to the degree of introduction of physical 
mechanism, the LIDAR RTMs can be divided into semi-empirical models, analyti- 
cal models, and computer simulation models (Table 2.1). The semi-empirical models 
use a digital surface model (DSM) or similar methods to represent the elevation dis- 
tribution of the Earth’s surface and simulates the LiDAR signal as a temporal sum- 
mation of Gaussian profiles that are computed by convolving the target reflectance 
and pulse energy distribution (Blair & Hofton, 1999). However, they only consider 
the single scattering of the laser on the target surface. Physically based analytical 
models use strong simplification of landscapes and ray propagations. For example, 
Sun and Ranson (2000) established a 3D LiDAR waveform simulation model, which 
simplifies the vegetation landscape into a series of 3D turbid voxels and only consid- 
ers the single scattering of laser rays. Ni-Meister et al. (2001) created a LIDAR wave- 
form simulation model based on 3D geometric optics radiative transfer (GORT), 
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TABLE 2.1 
LiDAR Radiative Transfer Model 
Category Model Proposer Features 
Semi-empirical DSM-based model (Blair & Hofton, 1999) DSM is needed as prior 
model knowledge; only single 
scattering is considered 
Analytical 3D LiDAR model (Sun & Ranson, 2000) Canopy is simplified as turbid 
model mediums; only single scattering 


is considered 

3D GORT model (Ni-Meister et al., 2001) Canopy is simplified as 
single-layer and double-layer 
forest; the ground is a flat or a 
terrain with fixed slope 


Computer RAPID (Huang & Wynne, 2013) Radiosity method 
simulation DART-RC (Gastellu-Etchegorry et al., 2016) Forward ray tracing 
model FLIGHT (North et al., 2010) Backward path tracing 

DIRSIG (Goodenough & Brown, 2017) Forward ray tracing 
LIBRAT (Disney et al., 2009) Backward path tracing 
RAYTRAN (Govaerts & Verstraete, 1998) Forward ray tracing 
FILES (Kobayashi & Iwabuchi, 2008) Forward ray tracing 
DART-Lux (Yang et al., 2022) Bidirectional path tracing, 


efficient and accurate 


which simplifies the ground as a terrain with a fixed slope or a flat surface without 
undulations and assumes that the ground has the same reflectance and the vegetation 
is a single-layer or double-layer forest. The semi-empirical and analytical models all 
make strong simplifications and assumptions on Earth’s scenery and laser propaga- 
tion process. Although computationally efficient, these models have an increasingly 
unsuitable accuracy for most potential LiDAR applications, such as signal modeling, 
parameter inversions of multiple targets, and LiDAR sensor design. 

With the development of computer graphics, some LiDAR computer simula- 
tion models based on Monte Carlo ray tracing (MCRT) (Pharr et al., 2016) were 
proposed, such as DART (Gastellu-Etchegorry et al., 2016; Yang et al., 2022) and 
FLIGHT (North et al., 2010). These MCRT-based LIDAR 3D RTMs are the most 
commonly used and accurate models since they are more adaptive to complex scene 
structures and usually do not make simplifications of ray interactions. According 
to the direction of ray tracing, they are usually categorized as forward models if 
rays are traced from the laser source, backward models if rays are traced from the 
receiver, and bidirectional path tracing (BDPT) models if rays are traced from the 
source and receiver simultaneously. Specifically, the forward ray tracing models use 
the laser source as the starting point for ray tracing and record the signal that finally 
enters the receiver. The backward tracing models trace the path from the receiver to 
ensure that the traced signal is received. The bidirectional tracing model traces the 
rays from the source and receiver synchronously. The purpose of all these ray tracing 
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models is to generate a series of random paths that connect the source and receiver. 
For LiDAR systems with the received field of view larger than the illumination area 
(footprint), the bidirectional tracing model usually has the highest accuracy and effi- 
ciency, followed by the forward tracing models and the backward tracing model with 
the lowest. 

Here we introduce a BDPT-based LiDAR modeling method called DART-Lux 
(Yang et al., 2022). The DART-Lux LiDAR model transforms the laser transport 
problem into an integration over all possible paths in the 3D landscape that con- 
nect the laser source and the receiver. It relies on the light transport equation, as in 
Equation (2.24). 


-dQ; 24 


L'u) = Ll Q)+ SLO, Q): flO" 2,2.) e086, 


Q 


where the exitant radiance L(r',Q,) at vertex r’ along direction Q, is the sum of 
the emitted radiance L, (r’,Q,) and the scattered radiance due to incident radiance 
L(r’,-Q;) along direction Q;, f (r’,-Q,; 505) is the bidirectional scattering distribu- 
tion function (BSDF) of the surface at vertex r’, and 0; is the incident angle between 
the incident direction Q; and the surface normal vector. Equation (2.24) can be trans- 
formed into an area integration form instead of a solid angle integration form, as in 
Equation (2.25). 


LO’ >= LE 3 r)+ pe >r) fr" 31 r) GO” >r’) ddr”) (2.25) 


A 


where r” and r are the previous vertex and next vertex of r’ in the light transport 
process, dA(r”) is the area at vertex r”, and the connection function G(r” © r’) 
between vertices r” and r’ is the product of |cos@,| and the Jacobian term to transfer 
solid angle integration over area integration. 

Since laser rays might be scattered multiple times in the 3D landscape, the received 
laser power can be represented by incrementally expanding Equation (2.25) to an 
infinite sum of a multidimensional integral, which can be expressed as a Lebesgue 
integral, as in Equation (2.26). 


Pumar (t(F)) = [ro du(r) (2.26) 
D 


where T(7) is the laser travel time along path 7 and ®,ipar (t(F)) is the returned 
laser power at the time T(7). D is the set of all paths. D=U%_,D, with D, the set 
of paths with n edges (n is also called the path depth). r is a path of the set of 
paths that connect the laser source and the receiver; it passes through a series of 
vertices in the scene: F €D = [717 =h, h,a, Fko, n EA, NE Zz}. f(T) is 
the power contribution of path r. du(7) is the area-product measure of path F, 
du(7,)=dA(r,)-dA(7,1)---dA(7% ). The solution of Equation (2.26) is the simulated 
LiDAR signal. The DART-Lux LiDAR model uses a Monte Carlo—based BDPT 
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FIGURE 2.5 Schematic diagram of the DART-Lux LiDAR model. (a) A sub-path starts 
from the LiDAR source and a sub-path starts from the receiver (at the same location as the 
source). (b) A series of paths that connect the LiDAR source and receiver are generated by 
the BDPT method. 


algorithm to solve Equation (2.26). Generally, the source and receiver of the LIDAR 
sensor are considered to have the same central position P, (x, y,H) and orientation 
Q, (8 LOL ). The paths from laser source to receiver are constructed based on the 
combination of two sub-paths that start from the source and the receiver, respec- 
tively, as in Figure 2.5. 

In the DART-Lux LiDAR model, the first random walk starts from the laser 
source and samples a light sub-path py, = po, P,---,Py,-1 With N, vertices, of 
which vertex po is on the laser source. Each vertex is sampled with the probabil- 
ity distribution function (PDF) pl pr) (k =0,1,...,N,-1). The second random walk 
starts from the receiver and samples a receiver sub-path gy, = qo.q.---.4y,-1 With 
N, vertices, of which vertex qo is on the receiver. Each vertex is sampled with the 
PDF plak ) (k=0,1,...,N,-1). The combination of the two sub-paths leads to the 
construction of a series of complete paths from source to receiver by using three 
sampling techniques: (1) “connect to light”: the vertex on the receiver sub-path is 
connected to a sampled vertex on the laser source; (2) “connect to receiver”: the 
vertex on the light sub-path is connected to a sampled vertex on the receiver; and 
(3) “connect vertices”: a vertex on the light sub-path is connected to a vertex on the 
receiver sub-path. Also, a light sub-path can randomly hit the LiDAR receiver, and a 
receiver sub-path can randomly hit the laser source. The length / (Ta) of any path s,, 
is the sum of the lengths of the light sub-path, the receiver sub-path, and the distance 
between the connected vertices p,_; and q,_;. The path length is further converted to 


the laser travel time T(r). The power contribution divided by the path PDF (2 3 j ) 


of a random path 7,, can unbiasedly estimate the laser returned power ripar. The 
simulated temporal power profile is the so-called LIDAR waveform. Furthermore, 
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the DART-Lux LiDAR model extends the single-pulse waveform LiDAR model- 
ing component to simulate multiplatform (satellite, airborne, terrestrial), multitype 
(waveform, discrete return, photon counting), and multipulse LiDAR signals. The 
multipulse simulation of laser scanning systems is an iterative loop of single-pulse 
modeling with specific geometry configurations and methods of detection and digi- 
tization. The point cloud of discrete return LiDAR is converted from the simulated 
waveform using Gaussian decomposition or peak detection methods. The points of 
photon counting LiDAR are derived per waveform using a statistical method based 
on the instrumental parameters of the single-photon detector. 


2.3. PRINCIPLES OF LiDAR ON DIFFERENT PLATFORMS 


According to different application requirements, LiDAR systems can be equipped 
on different platforms, including spaceborne (satellites and the International Space 
Station [ISS]), airborne (manned and unmanned aerial vehicles), ground (car, back- 
pack, terrestrial tripod), etc. This section introduces the principles of LiDAR on 
different platforms. 


2.3.1 Spaceporne LIDAR 


Until 2022, the existing spaceborne LIDAR systems for Earth observation include 
Ice, Cloud, and Land Elevation Satellite/Geoscience Laser Altimeter System (ICESat/ 
GLAS), Advanced Topographic Laser Altimeter System (ICESat-2/ATLAS), Global 
Ecosystem Dynamics Investigation (GEDI), and Gaofen (GF)-7. Among them, the 
ICESat-2/ATLAS implements the micro-pulse photon counting LiDAR technique, 
which is introduced in Section 2.4.3. The other systems use the full-waveform LiDAR. 
Here we mainly introduce the principles of spaceborne full-waveform LiDAR. 

The spaceborne LiDAR, with high platform altitude and repeated, large-area 
observations, can provide reliable and multitemporal height measurement data, 
which have been widely used in mapping the surface vertical structure of the Earth 
and the surface topography of other planets. For example, the first LiDAR satellite 
ICESat/GLAS system detects the Earth’s surface by emitting the laser beam, and 
each beam forms a nearly round footprint with a diameter of about 70 m on the 
ground (Schutz et al., 2005). The GEDI system, whose main mission is to monitor 
global forest biomass, carbon cycling, and biodiversity, has a footprint diameter of 
about 25 m on the ground (Dubayah et al., 2020). The GF-7 laser altimeter system 
has a footprint diameter of about 17 m. 

The footprint diameter of spaceborne LiDAR is usually very large (>10 m) due 
to the high observation platform. The large footprint diameter might cause the 
ground vertical extent to overlap with the object vertical extent. The degree of 
overlapping is affected by the terrain slope, sensor observation orientation, and 
footprint diameter. Usually, to avoid signal aliasing as much as possible, the space- 
borne LiDAR is designed with an on-nadir viewing mode, which means that the 
laser pulse is emitted vertically to the Earth’s surface. Even so, sometimes, in order 
to collect more altimeter data, some spaceborne LiDAR systems are designed with 
multiple laser ground tracks by beam splitting and dithering, e.g., GEDI and 
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FIGURE 2.6 Pulse geometry of spaceborne LIDAR for nadir and oblique orientations, with 
the illumination and reception areas at the horizontal plane. 


GF-7 systems. The modulation of multibeam splitting and dithering would cause 
the emitted laser pulse to be slightly incline with the milliradian level, which is 
usually negligible. 

Figure 2.6 shows the geometry configuration of the laser pulses for nadir and 
oblique views. Generally, the LiDAR transmitter (laser) and receiver are considered 
to have the same location P,(x,y,H) and the same center orientation Q L (01.91). 
0, and @, are the zenith angle and azimuth angle of the LIDAR orientation, respec- 
tively; H is the sensor altitude. The laser can be regarded as a point source that emits 
laser power into a fixed conical solid angle with the divergence half-angle of 0. The 
LiDAR receiver is usually a circular aperture with a radius of 7« and the field of view 
(FOV) half-angle of Orv. 

For the nadir view, the zenith and azimuth angles of the laser pulse are both 0°. 
A circular illumination area (i.e., Cim in Figure 2.6, also called the footprint) would 
be generated on the flat ground with an elevation of 0. The FOV area of the LIDAR 
receiver is also circular (i.e., Crp in Figure 2.6). The radii of the illumination area (rẹ) 
and FOV area (roy) are calculated as in Equation (2.27). 


lip = H: tan Of, Toy = H -tan Oroy + Nel (2.27) 


For the oblique view, the zenith angle and azimuth angle of the laser pulse are 
6, and @,, respectively. The illumination area (i.e., Eim in Figure 2.6) and the FOV 
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area (i.e., Erp in Figure 2.6) of the LIDAR system on flat ground are both elliptical. 
The long axis of the ellipse is in the same direction as the laser azimuth. In this case, 
the semi-major long axes (parallel to the laser azimuth) and the semi-minor axes 
(perpendicular to the laser azimuth) of the illumination area and the FOV area are 
calculated as in Equation (2.28). 


H - tan Oy H - tan 0y 
tib = a e ipen aa 
cos; cos“ 0, 
H - tan Oio (2.28) 
rh 1 
H - tan Op cos 0, E 
FL fov = Nel» N,fov = 
cos0, cos 0, 


where r,s, and 7,4 are the semi-major and semi-minor axes of the illumination area, 
respectively; 7, fov and 7,4, are the semi-major and semi-minor axes of the FOV area, 
respectively. 


2.3.2 AIRBORNE LIDAR 


Airborne LiDAR is defined as a LiDAR system equipped on the flying platforms 
such as drones and helicopters. It is also called airborne laser scanning (ALS), air- 
borne laser altimetry (ALA), airborne laser topographic/terrain mapping (ALTM), 
and airborne laser mapping (ALM), among others. Although the names are slightly 
different, the principles are similar. That is, the LIDAR senor generates laser pulses 
and receives returned signals, and the scanning device controls the orientation of 
the LIDAR sensor. The main features of an airborne LIDAR system are as follows: 


1. High precision. The vast majority of airborne LIDAR systems on the mar- 
ket use the pulsed laser ranging system with the consideration of its long 
measurement distance. However, it should be noted that the pulsed LIDAR 
requires advanced hardware and processing methods to achieve high rang- 
ing precision. Additionally, the ranging precision is influenced by the work- 
ing environment to some extent. The LIDAR usually works best in dry, 
cold, and clear atmospheric conditions. In contrast, the LIDAR performs 
poor in the daytime with strong sunlight. 

2. High power. The airborne LiDAR system scans the Earth’s surface in the 
air. It requires a relatively high laser power, so that the returned laser energy 
through long-distance atmospheric loss and target absorption is strong 
enough to be recorded by the detector. 

. Small size. The aircrafts carrying capacity and volume are limited. It is 
required to load LiDAR equipment and carry operators in a limited space. 
Hence, airborne LiDAR systems are commonly designed to be small and 
lightweight. 

4. Suitable wavelength. It would be better that the selected laser wavelength 

satisfies the following requirements: (1) the wavelength is located in the 
atmospheric window, so that the energy is less absorbed by the atmosphere; 


W 
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(2) the wavelength is positioned where target has strong reflectance, so that 
the target can return a strong laser signal; (3) the wavelength is positioned 
where the detector sensitivity is high, so as to reduce the entry of light in 
other bands; and (4) it is safe for human eyes. Most existing airborne LIDAR 
systems usually use a 1064-nm, 1550-nm, 905-nm, or 532-nm wavelength. 
Some use other special wavelengths, but relatively few. 


The airborne LiDAR systems scan the target in specific modes (Wehr & Lohr, 
1999). The four main scanning modes are as follows: 


1. Swing mirror scanning mode (Figure 2.7). This scanning mode drives the 


mirror with a motor to repeatedly swing within a certain angle, to achieve 
a series of scanning lines of laser beams on the ground. In one swing mir- 
ror cycle, the laser forms a periodic movement track on the ground, usually 
Z-shaped or sinusoidal. The advantage of the swing mirror mode is the 
simple principle, high scanning efficiency, and large scanning angle range. 
One of its significant disadvantages is that the laser spot density is not even, 
being denser at both ends and sparse in the middle, caused by the mechani- 
cal device continuously undergoing acceleration and deceleration processes 


(b) (c) 


FIGURE 2.7 Swing mirror scanning mode of airborne LiDAR. (a) Scanning diagram; (b) 
Z-shaped ground scan line; (c) sinusoidal ground scan line. 
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FIGURE 2.8 Rotating prism scanning mode of airborne LiDAR. (a) Scanning diagram; (b) 
parallel ground scan lines. 


during the swing periods. Also, this mode has the disadvantages of poor 
seismic performance and high requirements for motor performance. 

2. Rotate prism scanning mode (Figure 2.8). This mode causes the polygonal 
prism with the motor to rotate, causing the direction of the reflected laser 
beam to reciprocate within a certain range and scan on the ground. The 
scanning lines formed are usually parallel lines. Its advantage is that the 
obtained pulse density is uniform, as the angular velocity of the prism rota- 
tion is constant. The disadvantage is its relatively low scanning efficiency. 

3. Elliptical scanning mode (Figure 2.9). This mode forms an elliptical scan 
line on the ground. With the movement of the sensor, a series of overlapping 
elliptical scan lines is presented on the ground. Its advantages include the 
simple scanning principle, good seismic performance, high ranging accu- 
racy, and similar influence of inclination angle and atmospheric transfer on 
all laser points. However, the sampling points of the elliptical trajectory are 
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FIGURE 2.9 Elliptical scanning mode of airborne LIDAR. (a) Scanning diagram; (b) and 
(c) elliptical ground scan lines. 
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FIGURE 2.10 Fiber scanning mode of airborne LiDAR. (a) Scanning system; (b) ground 
scan lines. 


unevenly distributed, which brings difficulties in terms of data processing. 
Also, the high repetition rate of the spots on the ground cause relatively low 
scanning efficiency. 

4. Fiber scanning mode (Figure 2.10). Its principle is that the laser is trans- 
mitted through the fiber and enters the fiber group in turn; the terminals 
of the fiber group are arranged in a row; the fibers in different positions 
correspond to different fields of view; the lasers are emitted from the fiber 
arrays, in turn, to scan the target. The advantage is fast scanning speed and 
evenly distributed scanning points. The even points greatly shorten the data 
processing time. 


2.3.3 TERRESTRIAL LIDAR 


The terrestrial LIDAR systems include fixed LIDAR and mobile LIDAR such as 
vehicle-mounted and backpack-based scanning systems. This section introduces the 
working principle of the terrestrial LIDAR system based on fixed tripod. 

The terrestrial LIDAR fixes the 3D laser scanner on a tripod or other base sta- 
tion. The laser pulse is transmitted, then reflected by the object, and finally returns 
to the receiver. The distance between the scanner and object is calculated by the 
time delay of emitting and receiving the pulse. Also, the transverse scan angle 0 and 
longitudinal scan angle @ of each pulse are recorded. Generally, the operator sets 
the transverse scan angle range A@ and longitudinal scan angle range A@ to scan the 
scene in different directions. Also, the transverse scan angle resolution 60 and lon- 
gitudinal scan angle resolution 69 can be set to adjust the pulse density of terrestrial 
LiDAR. The number of emitted laser pulses (Nus) can be calculated according to the 
scanning angle range and resolution, as in Equation (2.29). 


A9 Ap 


Na =. 2.29 
s= 9 Bp (2.29) 
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Note the point cloud acquired by the terrestrial LIDAR is usually not uniform. 
Specifically, the objects close to the scanner have large point densities; otherwise, 
the point density is relatively small. Additionally, due to the occlusion between the 
objects, the single-station scanning cannot acquire comprehensive 3D information 
of all the objects in the scene. Therefore, in practice, the operators usually conduct 
multistation scanning with the terrestrial LiDAR and then splice the acquired data. 
More information about terrestrial laser scanning is provided in Section 3.2. 


2.4 PRINCIPLES OF LiDAR WITH DIFFERENT METHODS 
OF DETECTION AND DIGITIZATION 


Full waveform, discrete return, and photon counting are three common methods 
of detection and digitization for LIDAR systems (Mandlburger et al., 2019). Figure 
2.11 shows how the three methods record the emitted and received signals when 
the laser penetrates a complex forest canopy. Among them, the full waveform 
LiDAR (fwLiDAR) and discrete return LiDAR (drLiDAR) record the returned 
signal through an analog-to-digital convertor, while the photon counting LiDAR 
(pcLiDAR) directly records the returned laser photons. For the two modes of ana- 
log electronic detection and direct detection, the received signal is composed of a 
valid laser signal and invalid noise signal. The noise signal is mainly caused by the 
sensor system and radiation from the sun. The difference is that the analog elec- 
tronics detector converts the received laser power into the voltage signal, producing 
the returned intensity that varies over time. Also, the LIDAR with analog-to-digital 
record mode generally emits relatively high-power laser beams. The received valid 
laser signal has thousands of photons, which greatly exceeds the detector’s internal 
noise and solar noise, i.e., high SNR. It helps to accurately identify a valid signal and 
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FIGURE 2.11 Principles of LIDAR with different methods of detection and digitization. 
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then calculate the distance between the sensor and target. In contrast, the photon 
counting mode employs a micro-pulse laser that emits a very weak laser beam. A 
single-photon detector is used to detect the returned photons and record their return 
time. However, the solar noise is usually obvious compared with the weak laser 
energy. This causes the low SNR of pcLiDAR. 

We present the detailed principles of the fwLiDAR, drLiDAR, and pcLiDAR 
systems next. 


2.4.1 Futt-WaAverorm LIDAR 


Full-waveform LIDAR uses an analog-to-digital converter (ADC) to digitize the 
returned laser power time series into a waveform. The waveform can character- 
ize the complete vertical structure of the target, such as the gray shaded signal in 
Figure 2.11 (Harding et al., 2001). The pulse width of the returned waveform is 
influenced by the width of the emitted pulse and the vertical extent of the target 
within the footprint. For a target with a simple structure, such as a flat surface 
or building roofs, we would observe one returned laser peak. For a target with a 
complex structure, multiple return peaks are formed. For example, the laser hitting 
a building edge would be reflected by the building roof and the adjacent ground, 
thus forming two returned peaks. For objects with gaps, e.g., vegetation canopy, 
the laser pulse would irradiate the multiple objects in the transmitted direction, 
including leaves, stems, branches, and understory, forming a waveform signal with 
complex shapes. This process demonstrates the capability of LIDAR to capture 
vegetation’s vertical structure. 

Generally, the fwLiDAR system provides the digital recordings of emitted and 
received waveforms. The valid returned signal is calculated by detecting the first 
waveform location and the last waveform location exceeding the detected threshold. 
Then the ranging information is derived from the time delay between the center 
of the emitted pulse and the center of a valid received signal. The commonly used 
methods of calculating the pulse center include the centroid method, mean value 
method, and Gaussian fitting method (Tang et al., 2016). 

Figure 2.12 displays the development of fwLiDAR systems. The earliest wave- 
form sampling experiment for land terrain detection was the Airborne Oceanographic 
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FIGURE 2.12 Existing full-waveform LIDAR systems. 
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LiDAR (AOL) system developed by the National Aeronautics and Space Administration 
(NASA) (Krabill et al., 1984). Subsequently, scientists developed multiple types of 
airborne fwLiDAR systems by combining the full-waveform LiDAR sensor with the 
airborne scanning device. The typical ones are Airborne Topographic Mapper (ATM) 
(Krabill et al., 2002), Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) 
(Harding et al., 2000), Laser Vegetation Imaging Sensor (LVIS) (Blair et al., 1999), 
and Experimental Advanced Airborne Research LiDAR (EAARL) (Nayegandhi et al., 
2005). 

The main purpose of the SLICER and LVIS systems is to measure the vegetation 
canopy structure. They are designed with a large LiDAR footprint with a diameter 
larger than 10 m, of which obtained waveform sampling can reveal the forest vertical 
structure at a large scale. In contrast, the purpose of the ATM and EAARL systems 
is high-resolution topographic mapping with a small footprint (diameter <1 m). In the 
forest region, the small footprint cannot cover a tree crown completely. Hence, the 
acquired waveform cannot characterize the large-scale vegetation vertical structure. 

In particular, the EAARL system uses a 532-nm laser and transmits the returned 
signal to channels with different ADC ratios. Its purpose is to adapt the significant 
variance of laser intensity returned from different scatter surfaces, such as land, 
vegetation, surface of the water, and bottom of the water. Hence, the EAARL system 
is capable of acquiring bathymetric data and land topographic data. If the returned 
intensity is not scaled properly, the recorded waveform signal would be saturated or 
distorted. 

In addition to airborne platforms, full-waveform LiDAR is also equipped on 
spaceborne platforms to measure the Earth’s 3D surface structure on a global scale. 
However, the existing LiDAR satellites equip the laser altimetry, rather than the laser 
scanner. The ICESat/GLAS is the world’s first satellite-based, full-waveform laser 
altimeter for Earth observation. It was launched by NASA in 2003 and stopped work- 
ing in 2009. In 2018, NASA carried the spaceborne full-waveform laser altimetry 
GEDI onto the ISS for global forest monitoring. China’s GF-7, launched in 2019, and 
Terrestrial Ecosystem Carbon Inventory Satellite-1 (TECIS-1) launched in 2022, are 
both equipped with full-waveform laser altimeters. Their missions are global topog- 
raphy mapping and forest biomass monitoring, respectively. More details about the 
spaceborne LiDARs for Earth observation are introduced in Section 3.3 of this book. 

Except for Earth observation, full-waveform LiDAR systems were also applied 
to other planetary surveying, such as the Mars Orbiter Laser Altimeter (MOLA) 
for Mars surface mapping (Smith et al., 1999), Messenger Laser Altimeter (MLA) 
for Mercury surface detection (Cavanaugh et al., 2007), and Lunar Orbiter Laser 
Altimeter (LOLA) for lunar surface detection (Smith et al., 2010). 


2.4.2 Discrete RETURN LIDAR 


In contrast to fwLiDAR, discrete return LiDAR records only a few locations in the 
waveform time series output by the detector. Threshold detection, recording the position 
exceeding the threshold in the waveform time series, is the simplest method of record- 
ing discrete return points. With this method, only one or a limited number of objects 
irradiated by the laser are detected, which greatly reduces the recorded data volume. 
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A crucial system parameter of the drLiDAR system is the maximum number of 
returns that can be recorded per pulse. The earliest drLiDAR can only record one 
point per pulse. The drLiDAR systems developed subsequently gradually record 
two, three, and up to five points per pulse. Taking as an example of Figure 2.11, 
the drLiDAR detector records three discrete points from the waveform time series. 
In addition to recording the 3D coordinates (X, Y, Z) of each discrete point, the 
intensity (Z) of each discrete point is recorded by using the peak detection method 
and the Gaussian decomposition method. The peak detection method records the 
peak center of the waveform series as the point bin location and the peak amplitude 
as the point intensity. The Gaussian decomposition method performs the Gaussian 
decomposition algorithm on the waveform series. The center of the decomposed 
Gaussian component is recorded as the point bin location, and the amplitude is 
recorded as the point intensity. The time delay between the point bin location and 
the center location of the emitted pulse is used to calculate the distance between 
the target corresponding to each point and the sensor. Then combined with the 
global positioning system (GPS) and inertial measurement unit (IMU) data of the 
sensor, the 3D coordinates of each point are derived. Through multiple pulse laser 
scanning, massive discrete points with 3D geographic coordinates are generated, 
usually called a “point cloud”. 

The drLiDAR systems are usually designed with a small footprint and high pulse 
repetition frequency. Figure 2.13 shows an example of a discrete point cloud in a for- 
est area acquired by a small-footprint drLiDAR system, recording up to four returns 
per pulse. The position, intensity, and corresponding return number of each discrete 
point depend on the spatial structure and optical properties of the surface objects. 
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FIGURE 2.13 Multiple returns of discrete return LIDAR. 
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According to the LiDAR equation (Equation 2.20), the return laser intensity is related 
to the intercepted area and reflectance of the corresponding object. Once the object 
property makes it generate a strong enough laser intensity to exceed the detection 
threshold, the detector records it as one discrete point. The first discrete point corre- 
sponds to the location of the waveform signal that exceeds the intensity threshold for 
the first time, which might come from the top of the vegetation canopy, or the interior 
of the vegetation canopy, or the ground. The last discrete point corresponds to the 
location of the detected waveform signal for the last time, which might come from 
the interior of the vegetation canopy or from the ground. Sometimes, some ground 
points can be detected in the vegetation-covered region. This is because the laser 
pulse can penetrate gaps in the vegetation canopy to reach the ground. The frequency 
of detecting the ground under the canopy depends on many factors, such as canopy 
gap fraction, laser scanning angle, beam divergence angle, footprint diameter, pulse 
repetition frequency, and ground reflectance. 


2.4.3 PHOTON CountTiING LiDAR 


Different from the first two detection modes, the photon counting LiDAR records 
the arrival time of each laser photon to achieve the sampling of the target vertical 
structure (Figure 2.11). Through detection of multiple pulses, the accumulated pho- 
ton points potentially construct the spatial structure of the Earth’s surface (Howland 
et al., 2011). 

Full-waveform LiDAR and discrete return LiDAR systems usually emit high- 
intensity laser pulses, thus obtaining the laser signal with a high SNR. In contrast, 
photon counting systems emit very weak micro-pulses, with only a small number of 
photons returning per pulse. Hence, a single-photon detector is required to detect the 
returned photons. The very narrow pulse width of the emitted pulse (<1 ns), low jitter 
of the detector, and high resolution of the timing electronics make the pcLiDAR 
reach decimeter-level ranging accuracy. Note that the single photon detector of pcLi- 
DAR usually has some dark counts, i.e., system noise. The dark counts are quantified 
by the parameter “dark count rate,” which is the detector’s thermal noise count rate 
without any light source. The dark count rate is usually very low compared with the 
frequency of received photon counts. 

In addition to good detection sensitivity and low dark count noise, the single 
photon detector can detect multiple photons for one pulse. That is, a few photon 
points are recorded for each pulse. This characteristic is determined by a key detec- 
tor parameter, dead time, which is defined as the period needed by the detector to 
recover from a detection event, before another photon can be detected. The dead 
time of a single photon detector is very short, so that the pcLiDAR can distinguish 
between consecutive photons returned from closely spaced surfaces, e.g., the returns 
from vegetation canopy. This characteristic of pcLiDAR not only can detect multiple 
photons for one pulse, improving the photon density, but also record the acquired 
data in the form of discrete photon points. 

When the photon counting LiDAR system works in the daytime, the solar radia- 
tions reflected from the Earth’s surface and clouds are the main noise source. The 
solar noise rate is influenced by the surface reflectance, solar irradiance, atmospheric 
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environment, and sensor configuration, among other things. Generally, there is a fil- 
ter with a very narrow bandwidth located before the pcLiDAR detector. Its purpose 
is to prevent the electromagnetic wave other than the laser band from entering the 
detector. Also, pcLiDAR systems are typically designed to have a small FOV, con- 
fining the received signal within the FOV region. 

To summarize, the pcLiDAR has the characteristics of low laser power, small 
receiving FOV, and high pulse repetition frequency. Compared with the analog-to- 
digital LiDAR, it has strong potential in measuring fine 3D structures of the Earth’s 
surface at the large scale. Also, the output of pcLiDAR is the discrete photon points 
with specific 3D spatial coordinates. This makes the existing point cloud algorithms 
of data processing, analysis and visualization, and the software potentially suitable 
for pcLiDAR data. 


2.5 INFLUENCE MECHANISM OF THE SKY ON THE LiDAR SIGNAL 


The LiDAR is inevitably affected by the surrounding environment during the mea- 
surement process, especially spaceborne and airborne LiDAR. Figure 2.14 shows 
the radiative transfer process of sunlight energy (red arrows) and laser energy (green 
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FIGURE 2.14 Radiative transfer process of sunlight and laser among the atmosphere, 
Earth’s surface, and LIDAR sensors. 
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arrows) among the atmosphere, Earth’s surface, and sensors. It is seen that in the 
transmitted process, the laser pulse interacts with not only the Earth’s surface but 
also the atmosphere layer. Also, the sensor receives some noise signals from the 
sunlight, including the solar radiation scattered by the atmosphere and surface. This 
section presents the influential mechanisms of the atmosphere and sunlight on the 
LiDAR signal. 


2.5.1 EFFECT OF ATMOSPHERE 


The LiDAR systems usually use a laser with a short wavelength. The atmosphere 
will absorb and scatter the laser signal, causing the LiDAR sensor to be very sensi- 
tive to the atmosphere (Qiang et al., 2000). In particular, rain, dust, fog, and haze in 
the atmosphere interfere with the laser signal. The main influences of atmosphere on 
the LiDAR signal include (1) energy attenuation: the atmospheric gas, aerosol, dust, 
fog, and rain would absorb or scatter the laser signal, causing the attenuation of laser 
energy, and (2) laser refraction: the uneven spatial distribution of the atmosphere 
causes the laser to be refracted rather than directly fired. Additionally, atmospheric 
turbulence affects the fluctuation of laser energy distribution on the beam cross sec- 
tion and the beam expansion and drift. Atmospheric absorption in the laser transmis- 
sion path causes variations of the density gradient and refractive index of air, i.e., 
atmospheric thermal halo effect. This results in the nonlinear thermal distortion of 
the laser beam. The effects of atmospheric attenuation and atmospheric refraction in 
the LiDAR transmission process are introduced in detail next. 


2.5.1.1 Atmospheric Attenuation 


A laser has excellent monochromaticity. With the assumption that the atmosphere 
is uniform or stratified uniform, the laser energy through the atmosphere layer is 
expressed by the Beer-Lambert law: 


1(A,z) = (A, 0)e- 0 (2.30) 


where J(A,0) and /(A,z) are the initial laser intensity and the laser intensity passing 
through the atmosphere with a thickness of z, respectively, A is the laser wavelength, 
and o(A) is the atmospheric attenuation coefficient related to the wavelength. 

The atmospheric attenuation includes two parts: atmosphere absorption and 
atmosphere scattering. Therefore, o (4) can be expressed as Equation (2.31): 


0(A)=6n+Km+0.+K, (2.31) 


where O m» is the gas molecular scattering coefficient, Km is the gas molecular absorp- 
tion coefficient, o, is the aerosol scattering coefficient, and K, is the aerosol absorp- 
tion coefficient. 

The absorption coefficient of a gas molecule and aerosol for laser light can be 
calculated based on the transmission theory of monochromatic light. When the 
laser pulse is transmitted in the lower atmosphere below 20 km, the width of the 
atmosphere absorption spectrum is mainly determined by the pressure broadening 
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caused by gas molecular collision. The gas molecular absorption coefficient of the 
low atmosphere layer (K,,,) is expressed as Equation (2.32): 


S YL 
Aai T (V—vy) +¥n eas 
where v is light frequency, vo is the frequency of the center of the laser spectrum, S is 
the integrated intensity of the laser spectrum, 7, is a half-width of the Lorentz spec- 
tral line, and S and y; are both influenced by atmospheric temperature and pressure. 
In the upper atmosphere above 60 km, the gas molecular absorption coefficient (K mp) 
is calculated by the Doppler-broadened line-shaped function, as in Equation (2.33). 
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where Yp is the Doppler line half-width. In the middle atmosphere at altitudes of 
20-60 km, the collision broadening and Doppler broadening exist simultaneously. 
The gas molecular absorption coefficient (Kmy) is given by the Voigt line shape, as 
in Equation (2.34). 
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The atmospheric scattering effect includes Rayleigh scattering by atmosphere gas 
molecular and Mie scattering by spherical or quasi-spherical particles such as rain- 
drops, fog droplets, haze, and aerosols. 

Rayleigh scattering occurs when the laser wavelength is much larger than the 
particle size, and the gas molecular scattering coefficient Om is expressed as in 
Equation (2.35). 
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where n is the particle refractive index, N, is particle density in the air, and 6, is the 
depolarization factor. Based on experience, Om is expressed as: 


Om = 2.677 X10" Pvt /T (2.36) 


where P is pressure and T is temperature. 

Mie scattering occurs when the size of atmospheric particles is comparable to the 
laser wavelength. Generally speaking, the Mie scattering is the aerosol scattering. 
The total atmosphere attenuation coefficient can be approximately expressed as: 
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where Vm is visibility (unit: km), which is defined as the maximum distance at 
which the human eye can distinguish the target, and the coefficient b is related to 
the visibility. For normal visibility (6 km < Vn < 20 km), b = 1.3; For good visibility 
(Vn >20 km), b = 1.6; For poor visibility (Va <6 km), b = 0.585. Vie: 

In general, the scattering coefficients of atmospheric gas molecules and aerosols 
obey the negative exponential distribution law of altitude. That is, with the increase 
of atmospheric altitude, the scattering coefficients decrease rapidly. For example, the 
aerosol scattering coefficient above 5 km decreases more than an order of magnitude 
compared with that on the ground. Additionally, in actual situations, the aerosol scat- 
tering is the main interaction process in the low atmosphere layer near the ground, 
while the contributions of gas scattering and aerosol scattering are comparable in the 
middle and upper atmosphere layers. 


2.5.1.2 Atmospheric Refraction 

The atmospheric refraction effect is defined as the phenomenon that when the laser 
passes through the atmosphere, the light transmission path is bent and the traveled 
distance increases, caused by uneven atmosphere density or existence of an atmo- 
spheric refractive index gradient. Its effect on LIDAR measurement is to increase the 
measured error of target bearing and distance. 

The atmosphere density varies with the altitude, so the atmosphere refractive 
index for light at different altitudes is different. The refractive index n is related 
to the laser wavelength 7, atmosphere temperature T, humidity e, and pressure P, 
which is generally expressed as Equation (2.38): 


n=1+N(A,T,P,e) (2.38) 


where N is the refractive index modulus with the unit of order 10-6. For the standard 
atmosphere (P = 1 atm,T = 288.15 K,e = 0), the atmosphere refractive index from 
visible to near-infrared band is expressed as Equation (2.39): 


No(A) = 272.5794 + 1.5932A5° + 0.0155" (2.39) 


where No (A) is the standard atmospheric refractive index in the wavelength A. 
The atmospheric refractive index N(A) in any atmospheric condition can be cal- 
culated from the standard atmospheric refractive index No (A), as in Equation (2.40): 


N(A)= No(Ax (2.8434 10° x2 -o.1127% 5} (2.40) 


The atmosphere attenuation effect in the LiDAR radiative transfer process was 
introduced earlier. For practical application, other complex effects, such as atmo- 
spheric turbulence and thermal halo, should also be considered. The propagation 
of the LiDAR signal in the atmosphere directly affects the sensor performance. 
Therefore, it is necessary to fully understand the atmospheric propagation character- 
istics of the laser, find ways to avoid or reduce the atmospheric effects, and select the 
suitable working wavelength and working mode for the LiDAR system according to 
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the laser propagation mechanism in the atmosphere to finally achieve the high atmo- 
sphere transmittance of the laser beam. 


2.5.2 EFFECT OF SUNLIGHT 


The existence of sunlight makes the LIDAR RS radiative transfer process a RS con- 
figuration of “two light sources (laser and sunlight) + one sensor” (Figure 2.15). The 
laser source is regarded as a point light source emitting laser power onto a fixed 
conical solid angle. The emitted laser only lasts for a short time. The sun emits solar 
energy from an infinite distance to the Earth’s surface. Without the consideration of 
atmospheric scattering, the sunlight reaching the Earth’s surface can be regarded as 
continuous parallel beams with the same direction. Once the sunlight irradiates the 
surface within the LiDAR receiving FOV, the solar energy might be scattered by 
the surface and then received by the LiDAR sensor. The influence of sunlight on the 
LiDAR signal is different for different types of LiDAR systems. Next we introduce 
its influences separately. 

The large-footprint LIDAR systems (footprint diameter >10 m, usually mounted 
on spaceborne and airborne platforms) have a relatively large laser illumination 
region and receiving FOV region. The solar energy scattered into the LiDAR sensor 
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FIGURE 2.15 Radiative transfer process of dual light sources (laser and sunlight). 
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is relatively high. In contrast, the small-footprint LiDAR systems (footprint diameter 
<10 m, usually equipped on spaceborne and airborne platforms) have a small illu- 
mination region and receiving FOV region. Only very little solar energy enters the 
LiDAR receiver. Hence, the sunlight has less effect on the small-footprint LiDAR 
signal. 

For fwLiDAR, the influence of sunlight is shown as a relatively constant solar 
background noise. This is because the persistence of sunlight makes its radiation 
transfer process to adopt a dynamic equilibrium state, and the solar energy entering 
the LiDAR receiver almost remains constant. The solar noise of a LiDAR waveform 
can be easily removed by setting a background noise threshold. 

For drLiDAR, the solar noise is extremely weak due to the small footprint. Also, 
drLiDAR usually records discrete point clouds from waveform data by using the 
threshold detection method, peak detection method, or Gaussian decomposition 
method. The constant and very weak solar noise does not affect the calculated 3D 
coordinates of these points. 

For pcLiDAR, the detector treats solar photons in the filter band the same as 
laser photons. That is, the recorded photons are a mixture of solar photons, laser 
photons, and dark noise photons. Due to the extremely low laser energy of pcLiDAR, 
solar noise has a significant effect on the pcLiDAR data acquired in the daytime. 
Furthermore, due to the persistence of solar irradiation, the solar noise photons exist 
over the vertical extent of the entire scene, greatly increasing the density of the pho- 
ton points. Therefore, in order to obtain accurate surface structural information, the 
solar noise photons must be removed by postprocessing algorithms. The denoising 
algorithms for photon counting data are introduced in Section 4.3. 


2.6 SUMMARY 


This chapter first introduces the ranging principle of LiDAR systems, including 
pulsed LiDAR and phased LiDAR. Then, the LiDAR radiation principle is intro- 
duced based on the LiDAR equation. The commonly used LiDAR waveform model 
and LiDAR radiative transfer model are also presented. After that, we systematically 
introduce the principles of LiDAR systems with different platforms (spaceborne/ 
airborne/terrestrial) and different methods of detection and digitization (full wave- 
form/discrete return/photon counting). Finally, the influence mechanisms of the sky 
environment (atmosphere and sunlight) on the LiDAR signal are analyzed. 


EXERCISES 


1. Why do most LiDAR systems employ the pulse ranging method? 
2. What is the ranging principle of the pulsed LIDAR? How is it different from 
that of the phased LiDAR? 
. Please write the LiDAR equation and explain the meaning of each 
parameter. 
4. Please list three LIDAR 3D radiation transfer models and introduce their 
core ideas. 
5. Please state the commonly used scanning modes of airborne LiDAR sys- 
tems and their advantages and disadvantages. 


W 
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6. Please explain the differences of the LiDAR systems with three different 
methods of detection and digitization, i.e., full waveform LIDAR, discrete 
return LiDAR, and photon counting LiDAR. 

7. How are the sizes of the LiDAR footprint and field of view (FOV) calculated? 
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